Coupled Magnetic Excitations in Single Crystal PrBa 2 Cu 3 06.2 
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The dispersion of the low-energy magnetic excitations of the Pr sublattice in PrBa2Cu306.2 is 
determined by inelastic neutron scattering measurements on a single crystal. The dispersion, which 
shows the effect of interactions with the Cu spin- waves, is well described by a model of the coupled 
Cu-Pr magnetic system. This enables values for the principal exchange constants to be determined, 
which suggest that both Pr-Pr and Cu-Pr interactions are important in producing the anomalously 
high ordering temperature of the Pr sublattice. Measurements of the Cu optic spin wave mode show 
that the inter-layer Cu-Cu exchange is significantly lower than in YBa2Cu306.2- 

PACS numbers: 71.28.+d, 74.72. Jt, 75.10.Dg, 78.70.Nx 



The anomalous properties of non-superconducting 
PrBa2Cu 3 06+iE (PrBCO) jlj] remain an outstanding 
problem in the field of cuprate superconductivity. Its 
magnetic properties are strikingly different from those 
of other members of the RBCO family (where R is a 
rare earth, Y or La). Long range antifcrromagnetic order 
in both the Cu and Pr sublattices along with semicon- 
ducting resistivity persist over the entire range of oxygen 
doping < x < 1, and the magnetic transition in the 
Pr sublattice occurs at a temperature Tp r an order of 
magnitude higher than for other rare earths in RBCO 
It is now generally believed that these anomalous 
properties are caused by a hybridisation of the Pr 4/ 
and O 2p orbitals. In particular, an influential model of 
the electronic structure of PrBCO || explains the ab- 
sence of superconductivity by a localisation of the holes 
in hybridised Pr-0 bonds. Given the likelihood that the 
electronic and magnetic properties of PrBCO, including 
the occurrence of superconductivity both depend on 
this hybridisation, detailed information on the strengths 
and symmetries of the magnetic interactions is desirable. 

In this Letter, we describe inelastic neutron scatter- 
ing measurements of the magnetic excitations of PrBCO 
made for the first time on a single crystal. Using this 
data and a spin wave model of the coupled Cu-Pr sys- 
tem we have independently determined the magnitudes 
of the Pr-Pr, Cu-Pr and Cu-Cu interactions. We find 
that the Pr-Pr interaction is unexpectedly strong and 
would alone lead to Tp r being substantially higher than 
the rare-earth ordering temperature in other RBCO. We 
have also identified the effects of a pseudo-dipolar compo- 
nent to the Cu-Pr coupling on the magnetic excitations. 
Finally, we show that the inter-layer Cu-Cu exchange is 
reduced by a factor of two in PrBCO relative to YBCO. 

The experiments were performed on a crystal of 
PrBa 2 Cu 3 06.2 of mass 2g prepared by top seeding a 
flux. The crystal mosaic was ~ 1° and Tp r = 13 K. 
Measurements on the triple-axis spectrometers TAS 6 at 



Ris0 and IN8 at the ILL were made at energy transfers 
up to lOmeV and 70meV respectively. A horizontally- 
focussing analyser was employed to increase the mea- 
sured signal at the expense of some resolution in the 
scattering vector Q. The crystal was mounted in a 4 He 
cryostat and measurements were made in the (h, k, 0) or 
{h, h, I) scattering planes. 

Previous neutron scattering measurements on poly- 
crystalline samples of PrBCOe revealed two peaks cen- 
tred at 1.7 and 3.4 meV at T — 5 K, which shifted to lower 
energies as Tp r was approached ||. These peaks corre- 
spond to excitations of the Pr ions subject to the local 
exchange and crystal fields. To investigate the Q depen- 
dence of these excitations we performed energy scans at 
T = 1.8 K at points in reciprocal space primarily along 
the principal symmetry directions of the magnetic zone. 
Fig. 1 shows representative scans at Q = (0,0,2.2) and 
(0.75, 0.75, 0). Two peaks, labelled A and B, are observed 
at energies consistent with the polycrystalline data. Fig. 
1 also illustrates the variation in intensity of the excita- 
tions with the direction of Q. The intensity of peak A 
is largest when Q is in the ab plane, and is reduced for 
Q||c. Peak B exhibits the opposite trend. The width of 
both peaks (1-2 meV FWHM) is larger than the energy 
resolution of the spectrometer (0.4 meV FWHM). The 
implied intrinsic broadening is consistent with relaxation 
due to hybridisation as suggested earlier from the spec- 
trum of polycrystalline PrBCO §. The use of a single 
crystal, however, shows that the observed broadening is 
not solely due to dispersion. 

Single crystals also enable us to study the dispersion 
of the excitations, which is a central objective of this 
work. The insets to Fig. 1 show the low-energy data cor- 
rected for the non-magnetic background and elastic peak 
measured on a similar-sized crystal of YBCOg.2- The line 
shows a function comprising two damped harmonic oscil- 
lator response functions (including detailed balance fac- 
tor) fitted to the data. No variation in the energy of peak 
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B was apparent to within the experimental uncertainty. 
The energy dispersion of peak A was obtained by fitting 
a single response function to the data for E < 3meV. 
The results are shown in Fig. 2. A number of scans were 
also made at Q = (h, h, I) for I ^ 0, but no measurable 
dispersion along the c-axis was observed. 

We also studied the excitations of the Cu sublattice. 
The Cu spin waves in antiferromagnetic YBCO have been 
described by Tranquada et al. Due to the bilayer 
structure of YBCO the spectrum consists of two acous- 
tic and two optic modes, each pair being split into in- 
and out-of-plane modes by the anisotropy (a few meV 
in YBCO). The optic mode is of interest because its en- 
ergy depends on the exchange coupling between adjacent 
layers, which could be sensitive to the electronic struc- 
ture of the Pr ion. To measure the optic mode at the 
magnetic zone centre a number of constant energy scans 
were made in the (h, h, 0) direction through the point 
(5,5, 6.75), close to where the optic mode has maximum 
intensity. Scans at 45 and 55 meV are shown in Fig. 3(a). 
All scans at energies above ~50 meV have a peak centred 
at h = \, which we interpret as being due to the optic 
mode. The amplitude of this peak exhibits a stepwise 
increase with energy, as shown in Fig. 3(b). The op- 
tic mode gap, taken from the mid-point of the step, is 
53 ± 2 meV. The corresponding energy for YBCO6.2 es- 
timated in the same way from the available data fij] is 
70 ± 5 meV. 

We now consider a model to describe the dispersion of 
peak A. The model is based on the single-ion states of 
Pr 3+ in the tetragonal crystalline electric field (CEF). 
We obtained these states by refining a model for the 
CEF using the observed energy spectra of polycrystalline 
PrBCOe The refinement gave a doublet ground state 
with a first excited state a few meV higher in energy, con- 
sistent with previous analyses J5|,|| . To take into account 
the magnetic ordering at T < Tp r we included a molec- 
ular field. This field lifts the degeneracy of the ground- 
state doublet, hence we interpret peak A in Fig. 1 as the 
transition within the exchange-field-split doublet. Peak 
B is then the transition from the ground state to the crys- 
tal field singlet. We chose the magnitude and direction 
of the molecular field to reproduce the average energy of 
peak A and the observed direction of the moment || at 
about 45° to the c axis fiofl . The observed variation of 
the relative intensity of peaks A and B with the direction 
of Q is consistent with the cross-sections calculated from 
the CEF wavefunctions. 

The simplest model for the dispersion of peak A is then 
given by the Hamiltonian: 

H Pl = H CEF + ijp r ^(Ji - (J,)) • (J, - (J,)) (1) 

m 

where Jp r represents an isotropic exchange interaction 
between nearest-neighbour Pr ions (assumed to be > 



consistent with the observed antiferromagnetic ordering 
of the Pr moments), J; is the angular momentum of a 
Pr ion and Hqef is the single-ion crystal field Hamilto- 
nian including the molecular-field term. The summation 
in (0) is over all pairs of nearest-neighbouring Pr ions. 
We use the pseudo-boson approximation to calcu- 
late the dispersion from (|l]). This method gives a sim- 
ple relation between matrix elements of the calculated 
crystal field states and operators corresponding to de- 
viations of the angular momentum from the mean-field 
values. Here we only include the transition within the 
doublet ground state, as the exchange coupling does not 
give appreciable mixing with the singlet level. The result- 
ing dispersion of the doublet transition has two branches 
(due to the antiferromagnetic order) whose degeneracy 
is lifted by the non-zero x component of the exchange 
field. In practice, only one of the two modes is actually 
observed at a given Q owing to the variation in intensity 
between magnetic zones. The higher energy mode has a 
much larger cross-section near magnetic reciprocal lattice 
vectors that coincide with nuclear Bragg reflections (e.g. 
(1, 1,0)) whereas the lower energy mode has maximum 
intensity at the reciprocal lattice points corresponding to 
magnetic-only reflections (e.g. (|,|,0)). This intensity 
variation is in agreement with that observed experimen- 
tally, providing further confirmation for the identification 
of peak A with the intra-doublet transition, and for the 
assumption of an exchange field at an angle to the c axis. 

As a first approximation, we then attempt to fit this 
dispersion relation to the experimental data by varying 
Jp r . The best fit is given by Jp r = 0.029 meV, and is 
shown by the broken line in Fig. 2. Although the gen- 
eral form of the dispersion is reasonably well described by 
this model involving only the Pr sublattice, there is clear 
evidence of a discrepancy near the magnetic zone centre. 
This deviation suggests that we need to consider the ef- 
fect of interactions between the Pr excitations and the 
spin-wave excitations of the Cu sublattice. Owing to the 
very high Cu spin-wave velocity, these interactions are 
only important close to the magnetic zone centres at low 
energies. We model the Cu-Pr coupling with the same 
anisotropic term that has been found necessary to explain 
the observed non-collinear magnetic structure P Jl^Jl3]| . 

#Cu-Pr = ^ ^2 K <xp(ij)S a iJ(3j (2) 

where a, (3 label components {x,y,z}, K(ij) is a ten- 
sor describing the pseudo-dipolar interaction (consistent 
with the tetragonal symmetry of the crystal), and S de- 
notes the spin of a Cu ion. A similar term was also used 
to describe the excitations in Nd2Cu04 H]. 

For the interactions between Cu spins (-ffcu), we em- 
ploy the model used to describe the spin excitations of 
antiferromagentic YBCO || . This incorporates nearest- 
neighbour couplings in-plane J11 and between planes Jj_, 
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and an out-of-plane anisotropy. The optic mode gap at 
the zone centre is 2-i/jjj J±. For YBCO6.2 the value of 
J||, determined by Hayden et al 0, is 125±5meV. Here, 
we use a value of J11 = 127 ± lOmeV recently obtained on 
our crystal of PrBC06.2 by inelastic neutron scattering 
p6[ , consistent with Raman data pl|. Combined with 
our measurement of the optic gap, this value of J|| leads 
to J± = 5.5±0.9meV. We note the surprising result that 
the value of J± for PrBCOg.2 is a factor of 2 lower than 
that for YBC0 6 .2 (Jx = 11.5 ± 1.5 [§). 

Above Tp r the Cu ordering is similar to that observed 
in YBCC>6, however the presence of the interaction given 
by (||) leads to a modification of the Cu magnetic struc- 
ture at temperatures below Tp r . The magnetic moments 
on neighbouring Cu planes are observed to counter-rotate 
about the c axis by an angle <fi ||,|l7|]. This twist can 
be accounted for in our model by minimising the con- 
tribution to the ground state energy from the J± and 
pseudo-dipolar terms. The result is 

sin4> = 8{J x }{S}K xy /J ± (3) 

where J x is the x-component of the Pr angular momen- 
tum. The more general problem of verifying that the 
observed magnetic structure minimises the energy of the 
combined Cu-Pr system has been considered in [jl2| . 

Having modified -ffcu to take into account this ground 
state configuration, the excitation spectrum of the com- 
plete Hamiltonian H — Hq u + Hp r + i/cu-Pr is obtained 
from the bilinear terms in the spin-deviation operator 
expansion in the usual way. As a result of the pseudo- 
dipolar Cu-Pr interaction the spectrum is modified in 
two ways. First, the bilayer twist 4> leads to a gap in the 
in-plane acoustic Cu mode (the out-of-plane mode is also 
raised and continues to lie a few meV higher). Measure- 
ments below Tp r at higher energies indicate that the gap 
is ~ lOmeV. Secondly, the two Pr excitation branches are 
both significantly affected at the magnetic zone centre by 
coupling with the Cu modes Jl8| . As would be expected, 
the Pr spectrum is little affected away from the zone cen- 
tres because of the steepness of Cu spin- wave branch. 

Of the four independent components of the tensor K, 
the major contribution at the zone centre is from K xy . 
The other components produce only relatively small ef- 
fects, not varying very strongly with Q. Given that the 
Pr-Pr interaction can account for the dispersion away 
from the zone centre satisfactorily, we next fit the data 
with the dispersion of the full model by varying the pa- 
rameters Jp r and K xy . A good agreement is obtained, as 
shown by the solid line in Fig. 2 which corresponds to the 
values Jp r = 0.025 meV and K xy = 0.30 meV (the statis- 
tical errors in the fitted parameters are better than 10%), 
although near the nuclear zone centre the predicted dis- 
persion is so steep that the limited resolution of the spec- 
trometer will limit the extent to which the model can be 
tested. Further support for the fitted value of K xy comes 



from (|) which, together with (S) = 1/2 and (J x ) = 1.9 
from the single- ion model, gives <fi = 22°, in good agree- 
ment with the observed value of (f> = 20±5° ||. However, 
we should also mention two remaining shortcomings of 
the model. First, the fitted values for the exchange con- 
stants only produce a molecular field at the Pr site about 
half the size of that required to produce the observed en- 
ergy splitting. Secondly, the crystal field model predicts 
a Pr moment of 2 \xb compared to the measured value of 
1.2 [lb ■ These discrepancies may indicate that the elec- 
tronic state of the Pr ion is not completely accounted for 
by the crystal field model, and that other effects such as 
covalency need to be considered. 

Before concluding, we consider the magnitude of K xy 
and Jp r in comparison with other related measurements. 
Although no previous attempts have been made to deter- 
mine the Cu-R coupling in RBCO compounds, a similar 
pseudo-dipolar term seems to be present in Nd2Cu04, for 
which a value of K xy — 0.075 meV was found, a factor 
4 less than found here for PrBCO H]. As for the size 
of Jp r , we can judge this by making a simple estimate 
from mean-field theory of Tp r assuming Jp r to be the 
only interaction present. This requires a calculation of 
the single ion susceptibility using the same crystal field 
parameters as previously and predicts Tp r « 7K. This 
value is surprisingly large, suggesting that the high value 
of Tp r observed experimentally is due not only to a sig- 
nificant Cu-Pr interaction but also to large value of the 
Pr-Pr exchange. 

The central result of the work described in this paper 
is the set of exchange constants K xy , Jp r and Jx, deter- 
mined for PrBCO for the first time. These three con- 
stants are important because they reflect exchange path- 
ways that depend on Pr 4/ states, and so contain infor- 
mation on the underlying electronic structure of PrBCO. 
The Jp r , K xy are found to be surprisingly large, and 
this shows that both the Pr-Pr and Cu-Pr interactions 
are relatively strong in comparison to the corresponding 
interactions in other cuprates. Jx, on the other hand, 
is found to be a factor of two smaller in PrBCO than 
in YBCO, indicating that Pr has a significant influence 
on the interlayer exchange coupling. These results pro- 
vide insight into the anomalous electronic and magnetic 
properties of PrBCO, and suggest that the calculation of 
exchange interactions in PrBCO could be a very valuable 
test of models for its electronic structure. Similar exper- 
iments on an oxidised crystal of PrBCO, for which the 
effects of hybridisation could be even more pronounced, 
are planned for the near future. 

We acknowledge the financial support of the EPSRC, 
and the EU TMR Programme enabling the experiments 
to be performed at Ris0. 
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FIG. 1. Neutron inelastic scattering from PrBCO and 
YBCO at 1.8 K, for Q=(0, 0,2.2) and (0.75,0.75,0). The 
elastic peak from PrBCO is shown by the broken curve (re- 
duced by a factor of 50). Fits to the PrBCO data corrected 
with the YBCO as a non-magnetic background is shown in 
the insets. 
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FIG. 2. Dispersion of the low-energy Pr excitation. The 
measured data, obtained by fitting as shown in Fig. 1, are 
plotted with open symbols. The error bars are the standard 
deviations from the fitting procedure, but in some cases the 
points are the result of averaging several runs and the errors 
are then correspondingly reduced. The broken line is a fit to 
a model for the Pr sublattice only, while the solid line is a fit 
to the full model including a Cu-Pr interaction. (0.5, 0.5, 0) 
and (1,1,0) correspond to the centre of adjacent magnetic 
zones, so that the zone boundary is crossed at (0.75, 0.75, 0), 
but only the low energy mode has a measurable intensity in 
the first zone, and the high-energy mode in the second. 
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FIG. 3. (a) Measurement of the optic spin wave in PrBCO 
by neutron inelastic scattering. Constant energy scans at 45 
and 55 meV through (0.5, 0.5, 6.75) parallel to (h, h, 0) (in re- 
ciprocal lattice units (r.l.u)) are shown, (b) The amplitude of 
the observed peaks as a function of energy transfer. 
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